Translational Symmetry Breaking in the Superconducting State of the Cuprates: 
Analysis of the Quasiparticle Density of States 
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Motivated by recent STM experiments on Bi2Sr2CaCu20g+i [1~5], we study the effects of weak 
translational symmetry breaking on the quasiparticle spectrum of a d-wave superconductor. We 
develop a general formalism to discuss periodic charge order, as well as quasiparticle scattering off 
localized defects. We argue that the STM experiments in F^S^CaC^Os+i cannot be explained 
using a simple charge density wave order parameter, but are consistent with the presence of a periodic 
modulation in the electron hopping or pairing amplitude. We review the effects of randomness and 
pinning of the charge order and compare it to the impurity scattering of quasiparticles. We also 
discuss implications of weak translational symmetry breaking for ARPES experiments. 



I. INTRODUCTION 

Puzzling properties of the high-T c cuprates have of- 
ten been attributed to the existence of competing in- 
stabilities, and proximity (or even coexistence) of sev- 
eral ordered states. Possible instabilities that have 
been discussed in this context include charge density 
wave (CDW) order, non-two-sublattice spin density wave 
(SDW) order [6-9], spin Peierls order [10], and orbital 
magnetism [11,12]. Neutron scattering experiments on 
Lai.6- a; Ndo.4Sr :I ;Cu04 [13], La a; Sri_ a ;Cu04 [14,15], and 
La2Cu04+2; [16] demonstrated the coexistence of mag- 
netism and superconductivity, while recent experiments 
on strongly underdoped YBa2Cu.3O6.35 [17] have seen ev- 
idence of CDW order coexisting with superconductivity. 
Particularly striking in this context are recent STM ex- 
periments on Bi2Sr2CaCu20g+5 [1,2,5] which see spatial 
structure in the tunneling density of states with a period 
of four lattice constants. This structure was originally 
observed in the experiments in a magnetic field by J.E. 
Hoffman et al. [1], and later also seen in zero field by 
C. Howald et al. [2,5]. Modulo certain experimental sub- 
tleties, these experiments can be thought of as measure- 
ments of the spatial Fourier component (at the ordering 
wavevector Q) of the energy dependent local density of 
states (LDOS) pq(c}. 

In this paper, we demonstrate that the energy depen- 
dence of pc}(e) provides important information about the 
nature of charge ordering in these materials. It allows us 
to separate simple charge density wave order, that has 
only the Hartree-Fock potential modulation, from the 
more unusual charge orders that involve modulation of 
the electron kinetic energy (dimerization) or the pair- 
ing amplitude (anomalous dimerization). For example, 
when Q = (27r/4,0) (as observed in slightly overdoped 
Bi2Sr2CaCu208+a ), there is a change of sign in Pq(c) for 
energies around 40 meV when dimerizations are present, 
but not in the case of a simple CDW. When several of 
the simple distortions are present simultaneously, we can 
understand the resulting PQ(e) as (roughly) a superposi- 



tion of the corresponding simple cases, since the induced 
/9q(e) is approximately linear in the order parameter for 
small distortions. Such a superposition is necessary to 
understand the experiments of [2,5]. 

This superposition principle also applies when we have 
potential modulation at several wavevectors, and /3 q (e) 
can be analyzed separately for each wavevector q. This 
is necessary, for example, when we have randomness that 
pins the charge order, so that the single particle poten- 
tial is not a delta function in momentum space but has a 
narrow distribution centered at the ordering wavevector 
Q. This leads to a finite /? q (e) for a range of wavevectors 
around Q and, as we discuss below, taking a reasonable 
value of the CDW correlation length reproduces well the 
"weak dispersion" of the CDW peak observed in [5] . Our 
analysis can be extended to systems with no charge order 
but, instead, with localized defects, e.g. impurities. In 
this case we have a potential that is essentially momen- 
tum independent and we find strongly dispersing peaks 
in Pq(e) for a wide range of wavevectors. Such peaks 
have been observed in the STM experiments in [3,4] and 
discussed theoretically in [18,19]. We provide a quali- 
tative comparison of the STM spectra for systems with 
disordered CDW and impurity scattering. 

It is common to discuss spin density wave (SDW) order 
as the primary competitor to superconductivity in the 
underdoped cuprates [6,8,20-22]. An order parameter 
for non-two-sublattice magnetism is 



5(r) =0e lQ = r + 0*e- 4Qsr , 



(1) 



where the complex-valued vector (f> acquires an expecta- 
tion value in a state with broken spin symmetry. When 
the SDW order in (1) is collinear, it has an associated 
spin singlet order parameter that only breaks transla- 
tional symmetry and can be described as a generalized 
charge density wave [22] 



5p(r) = ipe lQ,:r + ip*e 



(2) 



Symmetry arguments determine the wavevector Q c = 
2Q S of such generalized CDW, but they do not clarify its 
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internal structure. For example, modulation of the local 
Hartree-Fock potential of the electrons and modulation 
of the electron kinetic energy (hopping) are both spin 
singlet order parameters that can be defined at the same 
wavevector and described by (2). Modulation of the elec- 
tron pairing amplitude also belongs to the same class of 
translational symmetry breaking since, in the supercon- 
ducting state with condensed Cooper pairs, order param- 
eters with charge two and zero are not orthogonal. It is 
important to note, however, that a long range SDW order 
is not a prerequisite for translational symmetry breaking. 
One can have a situation where quantum or thermal fluc- 
tuations destroy the spin order but preserve a long range 
order in the charge sector [23]. This was observed, for 
example, in underdoped YBa2Cu3 06.35 [17], where neu- 
tron scattering found period eight CDW but no static 
spin order. For slightly overdoped Bi2Sr2CaCu 2 08+<5, on 
which most of the tunneling experiments have been per- 
formed, neutron scattering experiments suggest dynamic 
spin fluctuations [24]. In our analysis we then assume 
that there is no SDW order and concentrate on the effects 
of spin singlet translational symmetry breaking. Another 
possible origin of a generalized CDW with no spin sym- 
metry breaking comes from pinning of SDW by disorder 
[19,25] or vortices [22,26]. 

We restrict our analysis to the case of weak transla- 
tional symmetry breaking, when the new order parame- 
ter can be treated as a small perturbation to the super- 
conducting mean-field Hamiltonian. This limit clearly 
applies to the experimental situation in [1-5], where the 
measured modulation is weak, and allows us to obtain 
explicit approximate expressions for /?Q(e). [This cir- 
cumvents solving a complicated set of equations numer- 
ically, as for instance carried out in [25].] Furthermore, 
we do not address the issue of the origin of charge order, 
but introduce it phenomenologically and study its conse- 
quences for the STM experiments. Our basic motivation 
is that a comparison of the energy dependence of /jq with 
experimental data can, in principle, be used to identify 
the correct order parameter(s) which, in turn, is crucial 
for understanding their origin. 

This paper is organized as follows. In Section II we 
introduce mean-field Hamiltonians that describe several 
kinds of translational symmetry breaking in a lattice sys- 
tem. For these generalized CDWs we derive an explicit 
formula for the Fourier component of the tunneling den- 
sity of states at the ordering wavevector. Section III dis- 
plays numerical results of this expression in the case of 
Bi 2 Sr2CaCu 2 08+ ( s type band structure and period four 
charge order. We show that recent STM experiments 
by [2,5] are consistent with a generalized CDW that has 
modulation of either the electron hopping or the pair- 
ing amplitude. We also consider period eight structure 
that may be relevant to YBa2Cu3 06.35- In Section IV 
we extend our analysis to phases with randomness and 
show that a realistic value of the CDW correlation length 



(20ao, with a>o the unit cell size) provides good agreement 
with the "weak dispersion" of the CDW peak observed 
in [5] . As a different application of our formalism we also 
consider localized perturbations in the crystal, such as 
impurity potentials, and demonstrate that these can ac- 
count for the strongly dispersing peaks observed in [3,4] 
at wavevectors not corresponding to the CDW order. In 
Section V we review how to include a more realistic model 
of the atomic wavefunctions, whose main effect is to in- 
troduce a momentum dependent structure factor. An 
important implication of this result is that the signal at 
wavevectors differing only by a reciprocal lattice vector 
should have peaks at the same energies, although their 
amplitudes may differ. We also discuss complications in 
the analysis of the STM data introduced by the normal- 
ization procedure used in the experiments. Finally, in 
Section VI we discuss several sum rules for the Fourier 
components of the density of states that may be useful 
for analyzing experiments. 

II. ORDER PARAMETERS AND MEAN-FIELD 
HAMILTONIANS FOR GENERALIZED CHARGE 
DENSITY WAVE PHASES 

Our starting point is a two dimensional one-band 
mean-field Hamiltonian that is commonly believed to be 
a good model for the physics of the d x 2_ y 2 superconduct- 
ing state of the cuprates: 

H = e k c L c k<x + A k (c kT c^ ki + c_ ki c kT ) (3) 

k<r k 

Here e k = — 2t(cos(k x )+cos(k y )) — 4ti cos(k x ) cos(k y ) — /i, 
A k = ^-(cos k x — cos ky) (from now on the unit cell size 
is set to unity), c rcr = A -1 / 2 ^ k c k(T e lkr , and A is the 
number of sites in the lattice. The Hamiltonian (3) can be 
diagonalized using the Bogoliubov transformation c k -\ = 
Mfc7feT + v kl-ki' c ~ k L = u kl-kl-v k -fl^ with u 2 k + vl = l, 
u k v k = A k /2E k , and E k = (e 2 + A 2 ) 1 / 2 . 

Weak charge order may be introduced into the state 
(3) by assuming the appearance of one or more of the 
additional order parameters listed below. We note that 
distinction between site and bond centered orders dis- 
cussed below is only defined for CDW with integer peri- 
ods. 

Site-centered charge density wave. The local Hartree- 
Fock potential is modulated along x with potential cx- 
trema on the lattice sites (see Fig. la): 

AWl = V k + C Lr C k+QoJ 

Bond-centered charge density wave. The local Hartree- 
Fock potential is modulated along x with the extrema of 
the modulation at midpoints of the horizontal bonds (see 
Fig. lb): 

AW 2 = Vo J2a, k ( a * c k+Q<T C k<r + a c k<T c k+Qff J , where 
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Longitudinal dimerization. Single electron tunneling am- 
plitudes are modulated on the horizontal bonds and the 
wavevector of modulation is along the same direction 
(i.e. along x). The bond centered version, in which 
the extrema of the modulation lie on the bonds (see 
Fig. lc) corresponds to: AHz = Vq E ct k cos (^a; + 



where 



k+Q(T c k(T + acj C(J c k+Q(T 



a.) 

2 1 

Transverse dimerization. Single electron hopping is mod- 
ulated on the vertical bonds, and the wavevector of mod- 
ulation is along the horizontal direction (i.e. along x). 
The site centered version (i.e. with extrema of the mod- 
ulation realized on the vertical bonds) is shown in Fig. 
Id, and corresponds to: 



AH 4 



VbEk<r cosk v 



k+Qo- c kcr + c k(T c k+Qo^ 

Anomalous longitudinal dimerization. The x-components 
of the d x 2_ y 2-w&ve pairing amplitudes are modulated 
in the x direction. The bond centered version, shown 
in Fig. lc, corresponds to: AH5 ~ Vo E<t k {cos(fc 2: + 

%) ( a * C ik+QT C -kl + a C kT C -k-Ql) + h.C.}. 

Anomalous transverse dimerization. The y-components 
of the d x 2_ y 2--w&ve pairing amplitudes are modu- 
lated in the x direction. The site centered ver- 
sion, shown in Fig. Id, corresponds to: AH$ = 

V o E CT , k { cos k v ( c k+Q T c -ki + c k T c -k-Qi) + h - c -}- 

Note that these subdominant order parameters may 
appear either as a result of a phase transition in the 
bulk, or due to pinning by vortices, impurities or any 
other defects (see discussion in Section IV). Following ex- 
perimental observations in [2,5,27], we assume that the 
order is unidirectional, and choose the ordering wavevec- 
tor to be Q = Qe x [28]. However, even if we were to 
assume checkerboard order, our analysis is carried out to 
linear order in perturbation theory and, by linear super- 
position, our results would be identical to those obtained 
for unidirectional order. For the STM experiments in 
Bi 2 Sr 2 CaCu208+5 [1,2], Q — 27r/4, while the neutron 
scattering experiments on YBa2Cu306. 35 [17] correspond 
to the smaller value: Q — 2ir/8. We point out that the 
six cases listed above are, in general, not orthogonal to 
each other in a symmetry sense. As a result it is con- 
ceivable that more than one order parameter could be 
simultaneously non-zero; for example, in a microscopic 
model without particle-hole symmetry, a simple CDW 
would be expected to induce dimerization as the two or- 
der parameters are linearly coupled [29] . 

Upon expressing the Hamiltonians AH in the basis 
of Bogoliubov quasiparticles, they reduce to the generic 
form 

Am = ]T [A l fe7 L7k+Qa + A% 7 k+QCT 7ka 

ka 

+ E [ B UU+Ql + B k*7 k+ Q T 7ki + h-c] , (4) 



A 1 

A 2 - 
^k - 



= V u) k 5 k = Vb?7 k 

Voaw k -B k = Voafy k 

Al = V a a cos(k x + §) uj k Bl = V a cos(k x + §) r? k 



At = Vo cos(ky) w k = V cos(fc y ) r? k 



(5) 



A k = -V a cos(k x 



§)?7k Bl = V a cos(k x 



A^ = -V cos(fc y ) rj k B^ = V cos(k y ) tu k 



-)^k 



in terms of the coherence factors w k = u k +QW k — vw+qv k 
and ?7 k = ?i k+Q « k + v k+ QW k . 

STM experiments measure the local density of states 
P(r,e) ^ = EnA\( n \4«\0)\ 2 5(e-e na ) + 

|(n|c r(T |0)| 2 (5(e + e„o)}, where the summation over n 
ranges over all excited states. In particular, we are inter- 
ested in the Fourier transform 



^ £ ) = ^E e " qr ^ r ' e ) 

[(0|c k+qCT |n) (n|4 CT |0)£(e - e n0 ) 
+ (0|4 CT |n) (n|c k+q(T |0)<5(e + e n0 ) 



N 

r 

-Y 

nkcT 



(6) 



Although a full treatment of all terms in (4) is compli- 
cated, progress can be made if we assume that the or- 
dering represented by AH is weak, allowing us to obtain 
an analytic expression for the Fourier transform that is 
exact to linear order in Vo- This is then the sum of two 
contributions 



P^)^P^(e)+pl(e)+0(V 2 



;> 



(7) 



where Pq(e) is obtained by ignoring the i3 k term in the 
perturbation (4) and vice versa. 

A small value of Vo leads to another important sim- 
plification: We only need to consider the pairwise mix- 
ing between states connected by AH. For instance, 
in computing /5q(e) for Q = 2n/A one would have 
to analyze coupled equations for four quasiparticles 
(k, k + Q, k + 2Q, and k + 3Q = k - Q) con- 
nected by the perturbation. However, in the limit when 
V) is small, there is at most one pair of quasiparti- 
cles that have similar energies, and that will be hy- 
bridized appreciably by AH. This hybridization can 
be analyzed by diagonalizing the corresponding two- 
by-two Hamiltonian, which gives the new eigenstates 
|a kCT ) and |/3 k(T ) with energies I? k ± = ' Ek+ f k+Q ± 

nl/2 



ik+Q \2 



) 2 + l^k 



Note that these states sat- 



isfy (0|c k+Q(T |a: k(J )(a! k(T |c k(T |0) = I u k u k+Q sin 26> k e 
(0|c k+Q(T |/3 k(T )(/3 k(T |cj c(j |0) = -iu k u k+Q sin26> k e- 
where we have defined A k = |^4 k |e lXk , and tan2# k = 
2|yl k |/(.E k — -E k +Q). From these results one easily finds 



A* 
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UkUk+q (<J(e-£k+) -5(e--Bk-)) (8) 
+Wk«k+Q (<S(e + #k+) - 5(e + i?k- 



When considering p||, one would naively expect that 
it is always smaller than Pq, because the perturbation 

terms of the form 7k±Q|7-ki connec t states that differ 
in energy by E^ + E^+q, a factor that is never small. 
However, in some cases the coherence factors in A k van- 
ish at important regions of the Brillouin zone, making 
Pq(e) anomalously small. In addition, as we discuss be- 
low, both p A ' B (e) are large at biases corresponding to 
the saddle points on the degeneracy lines E^ = i?k+Q 
and van Hove singularities of the Bogoliubov quasiparti- 
cles e w A . A nearly identical analysis of the one above 
for p A yields 



N ^ 

k 



+ \Bk\ 



WkWk+Q<5(e - E k+ ) + u k +QWk<5(e - #k-) 
-Uk+QVu_5(e + E k+ ) - M k Wk+Q^(e + £ k _) 



where E^± = ± 



-Ek — -Ek+Q 



Ek + Ek+Q\2 



(9) 
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Equations (8,9), are two key results of this paper. In com- 
bination with Eqns (5), they provide an explicit expres- 
sion for the energy dependence of the Fourier component 
of the local density of states pq (e) when the translational 
symmetry breaking is weak. 

From the form of Ak and Bk, it is obvious that when 
there is no mixing between bond and site centered CDW, 
Pq (e) can be made real at all energies by an appropriate 
choice of the overall phase, i.e. by a shift in the origin of 
coordinates when doing the Fourier transform. One ob- 
vious observation is that the results for the site-centered 
and bond-centered CDW are identical modulo an overall 
phase factor of e 1 ®/ 2 . If one defines the Fourier trans- 
form in such a way that it is real in both cases, the origin 
will coincide with one of the sites of the lattice for the 
site-centered CDW, and it will be at the center of a bond 
for the bond-centered CDW. Hence careful analysis of the 
STM data allows one to distinguish two kinds of CDW, 
a task that is not possible in neutron scattering exper- 
iments with current resolution. Mixing site and bond- 
centered orders breaks inversion symmetry and leads to 
a complex- valued pcj(e). 

III. CHARGE ORDER WITH NO RANDOMNESS 

A. Period four CDW in BiaSraCaCuaOs+s 

We first focus on modulations at Q = (27r/4,0) that 
is relevant to Bi2Sr2CaCu208+5 [1,2]. Figure 2 shows 



results of the numerical evaluation of formulas (7), (8) 
and (9) for various perturbations (5). [As transverse and 
longitudinal dimerization curves are qualitatively similar, 
curves corresponding to the former are not displayed.] 
We choose the band structure and the value of A in (3) 
appropriate to Bi 2 Sr2CaCu208+(5 : h/t = —0.3, p/i = 
-0.99 (this corresponds to 14% doping), A /t = 0.14 
and A = 40 mcV [31]. We set V /t = 0.02, although its 
precise value is inconsequential, as pcj(e) scales linearly 
with Vo when the latter is sufficiently small. 

If we turn our attention to the expression for p A , 
Eqn. (8), we see that the energy denominator is small- 
est for those quasiparticles lying close to the degeneracy 
points i?k = iik+Q, which are strongly hybridized by the 
Ak part of the perturbation. Figure (3) shows the four 
loci of such points, a through d, that are degenerate with 
a' to d! respectively. The pairs aa! and W are obvious, 
since they have k x = ±7r/4 and k x — ±37r/4 (for the 
same k y ); the other two require a more detailed analy- 
sis of the band structure. Out of the set of degeneracy 
points, we expect large contributions from the neighbor- 
hood of points A and B, as the dispersion of hybridized 
energies Ek± is fiat at these points. These same regions 
of the Brillouin zone will dominate the ps contribution, 
since the energy denominator in (9) will be small only if 
both k and k + Q lie close to the Fermi surface, which 
occurs only in the neighborhood of points A and B. In 
addition we expect, for both p A and p B pieces, a large 
contribution at e = Ao, where a van Hove singularity for 
the Bogoliubov quasiparticles yields a logarithmic diver- 
gence in the density states. 

We turn now to the numerical results displayed in 
Fig 2. Consider first the simple CDW curves. The sharp 
features that dominate the CDW plots can be understood 
in terms of the degeneracies mentioned above: the peak 
at energies around 0.5Ao comes from the vicinity of the 
A point, the peak around 0.7Ao comes from the vicinity 
of B, and the pile around Ao comes from the van Hove 
singularity near the (0, it) and (7r, 0) points. The longi- 
tudinal dimerization results can be similarly understood 
by taking into account the additional minus sign in the 
vicinity of the point B due to the cos(k x + 7r/4) factor 
in Ak and -Bk- The results for the anomalous dimer- 
ization can also be understood in this framework after 
taking into account the extra sign modulation in itk«k, 
which changes sign whenever Ak does. Note that, for 
all perturbations considered, Pq(c) displays approximate 
particle-hole symmetry for small biases, as observed in 
STM measurements [2] . This is not a generic property of 
pcj(e); for example, for a diamond-shaped Fermi surface 
p = t' = the CDW signal is exactly antisymmetric. 
Finally, note that pq(e) goes to zero at e = in all cases; 
this reflects the vanishing density of low-energy quasipar- 
ticle states in an ideal ci-wave superconductor. 

While the results in Figure 2 describe a system with 
infinite quasiparticle lifetime and no disorder, in a real 
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system disorder will smear the sharp features in Pq(c). 
To model this, these curves are re-displayed in Figure 4 
after smearing over an energy width w — 0.2A . This 
procedure smooths the sharp features in the spectra, and 
generates finite intensity at low energies. Notice that the 
smeared CDW curve does not have the two large peaks 
surrounding zero bias, nor does it have clear zero cross- 
ings at |e| « A , the dominant features of the STM spec- 
tra observed in [1,2]. By contrast, the signal related to 
longitudinal dimerization or, especially, to either kind of 
anomalous dimerization, share many of the qualitative 
properties of the data. However, neither curve by it- 
self accounts for all the observed features in the data. 
This prompts us to consider a combination of several 
kinds of order. For example, if we assume that no pair- 
ing modulation is present, the combination of longitu- 
dinal dimerization and CDW, (long.dim.) + 1.05(CDW), 
shown as a solid curve in Figure 5 reproduces the STM 
results reasonably well, with only a small difference in 
the position of the peaks (±17 meV, compared to ex- 
perimentally observed ±25 meV). Alternatively we can 
match experimental data by considering the combina- 
tion of anomalous longitudinal dimerization and CDW, 
(anom. long. dim) + 0.2(CDW), shown as a dashed curve 
in Figure 5. It slightly overestimates the peak bias to be 
±29 meV, and yields a low intensity at zero bias. Any in- 
termediate combination between these two scenarios also 
gives good agreement with experiments. Although CDW 
was used in both combinations discussed above, it can 
be substituted by transverse dimerization, which yields 
a qualitatively similar pq(e) to CDW. We note that, for 
e < 3Ao, the results come from the vicinity of the Fermi 
surface and are robust against variations in the band 
structure that do not alter qualitatively the shape of the 
Fermi surface (e.g. the a and b' lines do not move below 
the Fermi surface). 

We note, however, that a certain care should be ex- 
ercised when comparing our results to the STM spectra 
in [1-5]. An additional complication of the experiments 
is that for every point on the surface of the sample the 
height of the STM tip is adjusted to keep the tunneling 
current at a certain voltage fixed. This implies that the 
local density of states is not measured directly, but in- 
stead its product with some space dependent function is 
measured. In Section V we review how this normalization 
procedure can be included in analysis. 

B. Period eight CDW in YBa 2 Cu 3 06.35 

To model YBa 2 Cu3 6.35 for which CDW-type peaks 
have been observed at Q = 2n/8 [17] we take the same 
band structure t\/t — —0.3, but a different value of the 
chemical potential fi/t — —0.815 (this corresponds to 
6% doping). We use the same value of Ao/t = 0.14, 
A = 40 meV, Vo/t = 0.02, and keep the energy smearing 



w = 0.2Ao. The main difference with the charge order 
at Q = 2-k/A is that the analog of line a in this case is 
inside the Fermi surface, so that the only contributions 
will come from the vicinity of point B at energies around 
O.8A0. This leads to less structure in pq(c) and smaller 
intensity at zero energy (see Fig. 6). 

IV. DISPERSION OF THE STM SPECTRA 

Recent experiments [3,4] demonstrated that the STM 
spectra of Bi 2 Sr2CaCu208+,5 cannot be explained by 
charge order at a unique wavevector. Peaks in p q (e) have 
been observed away from (27r/4, 0) and the wavevectors of 
the peaks are energy dependent. In this section we review 
and compare two possible scenarios for such dispersion of 
the STM spectrum: 1) randomness and pinning of charge 
order, 2) scattering of BCS quasiparticles by impurities 
and crystal defects. Both cases can be described using 
an extension of the formalism presented in the previous 
section. We consider a single particle Hamiltonian that 
generalizes equation (4) 

AW = ^[y kq 4 >ff c k+qi(7 + h.c. (10) 

kqcr 

+Wkq (c k ,<x ct -k-q,- CT + C_ k ,_ ff C k+qi(7 ) + h.C. . 

Here q describes the wavevector of the potential modula- 
tion, and the k dependence of V and W gives its internal 
structure (e. g. simple CDW vs dimerization) [32]. In Sec- 
tion II we considered charge order at a unique wavevec- 
tor that corresponds to taking potentials V and W as 
S(q — Q). In the case of a disordered CDW we expect 
that these functions are no longer ^-functions but are cen- 
tered narrowly around some particular wavevector. By 
contrast, when translational symmetry breaking comes 
from impurities, we expect to find V and W that extend 
over a wide range of wavevectors q. A crucial property of 
(8) and (9) is that p q (e) depends linearly on the pertur- 
bations Vq and Wq, hence the formalism for computing 
p q (e) can be applied independently to each wavevector 

q- 

The charge order observed in Bi 2 Sr 2 CaCu208 + 5 [2,5] 
had strong signatures of randomness and pinning in 
the form of lattice defects. The correlation length 
estimated from the distance between defects was ~ 
20ao. If we assume the charge order to be of the form 
1.05(CDW)±(long.dim.), we can describe it as 

Kkq = Vo(q)(l.05 + cos(A: x + |)), IT kq = (11) 

where Vo(q) is a Gaussian distribution function centered 
at (27r/4,0) with a width 27r/20a . We display in Fig. 7 
the signal produced by a perturbation of this kind, for 
bias voltages 8, 12, 16, and 20 mV, as a function of 



5 



wavevectors along the (0,0) to (ir, 0) direction. The re- 
sulting dispersion agrees closely with that observed in [3] 
and [5]. 

In experiments of Hoffman et. al [3] and McElroy et. al 
[4] peaks in the LDOS were observed at very different 
wavevectors from (27r/4, 0) (including some in diagonal 
directions). This suggest that either Vk q or Wk q must 
be non-zero over a fairly wide range of values of q, and 
the most natural candidate is scattering by impurities 
[3,4,18]. For concreteness, we assume that the impurity 
induces a higher chemical potential at a single site, so the 
perturbation used corresponds to a simple CDW which is 
uniform in q, Vk q = Vq, Hkq = 0. In Fig. 8 we show the 
signal computed along the (0, 0) to (ir, 0) direction at bias 
voltages 8, 12, 16, and 20 mV. In all cases there is a pro- 
nounced peak that disperses with the applied bias volt- 
age. To find the positions of these peaks we reverse the 
arguments given in Section II. There, we started with a 
potential at wavevector Q and found that only quasipar- 
ticles at certain energies were strongly affected by it. Now 
we need to find the modulation wavevector that affects 
quasiparticles at a given energy. From the band struc- 
ture of I^S^CaC^Os+a we find for the peak positions 
(in units of 2tt): 0.35, 0.32, 0.29, and 0.26. The curves 
on Fig. 8 show general agreement with this "quasipar- 
ticle scattering" argument [3,4], except for a consistent 
small shift to lower wavevector, which comes from the 
energy smearing procedure. This dispersion is stronger 
than that displayed in the data at wavevector (27r/4,0), 
but is in good agreement with the dispersion observed at 
other wavevectors. 

In the discussion above we considered two situations: 
ordered CDW and non-interacting electrons with impu- 
rities. There may also be an intermediate regime with 
interacting electrons close to the CDW instability and 
with disorder [19,30]. Qualitatively, this case may be de- 
scribed by equation (10) but with the potentials Vk q and 
Wk q coming not only from the external fields but also 
from the density induced in the electron system. For 
simplicity let's take only one of the channels discussed in 
Section II, e.g. simple CDW (the generalization to the 
case of several channels is straightforward). Then 

AH = £ p q (7 q/ 5 q + J2 K Xt Pl + h - c - ( 12 ) 
q q 

with 

Pi = X! C ko- C k+qcr 

Response of the quasiparticles to (12) in the Hartree ap- 
proximation is determined by the effective perturbation 
Hamiltonian 

AH = ]T V^ +qff c k(T + h.c. 

kqtr 



Vf = + t/ q (p q ) 



T/CXt 

*q 



i - t / qXo(q,w = 0) 



(13) 



where xoiQi -' = 0) = (/5 q /5 q ) must be computed for the 
non-interacting system. This corresponds to contribu- 
tions to the LDOS from the class of diagrams in figure 
9(a). When the system is close to the CDW instability, 
the denominator of (13) approaches zero around some 
particular wavevector Q. Hence, V^ s may be peaked 
around Q even when Vq Xt is momentum independent. In 
principle we can go beyond the Hartree approximation, 
including diagrams such as those shown in figure 9(b,c). 
These will introduce a frequency-dependent self energy, 
as in the case of pinned spin density fluctuations consid- 
ered in [25]. 

It is interesting to ask whether by analyzing experi- 
mental data one can separate contributions of disordered 
CDWs from those due to quasiparticle scattering off im- 
purities. Reference [5] pointed out that the presence of 
a weakly dispersing signal at wavevector (27r/4, 0) makes 
charge order a likely candidate at that wavevector. Here 
we build on this idea and suggest that a more consistent 
approach is to analyze /o q (e) at many different wavevec- 
tors using a reasonable set of basis functions, e.g. simple 
CDW and dimerization. Such analysis will give the q 
dependence of various components in the potentials Vk q 
or Wk q - We expect that some of them will be almost 
uniform in momentum space and correspond to localized 
defects, such as impurities; whereas others will be cen- 
tered around particular wavevectors and arise from the 
existence or at least proximity of charge order. 

When the field of view of the STM measurement con- 
tains more than one impurity there are several impor- 
tant questions that we need to address. We must ask 
how the contributions from different impurities add up 
and whether the system retains the property of unifor- 
mity of phase of p q (e) at fixed q but different e. We 
consider impurities that cause an arbitrary potential Vkq 
in Eqn. (10), i.e. they may modify the chemical po- 
tential, the electron kinetic energy, or the pairing ampli- 
tude, but first we assume that all impurities are identical. 
The Fourier component of the LDOS is proportional to 
V^°* = Vkq J2r a e * qra ! wnere r a runs over impurity posi- 
tions. If one impurity does not break parity symmetry, 
we can make Vk q real by choosing the origin at the po- 
sition of this impurity (see also discussion in Section II). 
This implies that V^°' has a phase that depends on q 
only and is the same for all k, which in turn proves that, 
at fixed q, p q (e) has constant phase (modulo ir) for all 
values of e (sec Eqns (5)-(9)). So, in the case of identical 
impurities we have only one phase to worry about and 
we can always make p q (e) real by an appropriate choice 
of origin. When impurities are different, we will have 
an intrinsically complex Pq(e), with possibly energy de- 
pendent phase at different bias voltages. In either case, 
interference among the impurities leads to an apprecia- 
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ble suppression of the amplitude of p q (e). When there 
are many impurities in the area A of the STM field of 
view, and their positions are uncorrelated, each impurity 
introduces a random phase to V^°* , whose amplitude can 
be analyzed in terms of a random walk. Therefore, in 
a typical experiment we expect that with increasing the 
system size, (|p q (e)|) will decay as l/^/A, with statisti- 
cal fluctuations of the same order. This argument also 
applies to the case of disordered CDW, where the role of 
impurities is played by defects in the CDW lattice. 



V. EXPERIMENTAL CONSIDERATIONS 

Our discussion in the earlier sections was restricted 
to models on a square lattice for which we calculated 
the lattice density of states in the presence of several 
kinds of translational symmetry breaking. In analyzing 
actual experimental data, additional effects need to be 
taken into account: the real space structure of the atomic 
wavefunctions, and the current normalization condition 
used in the STM measurements [1-5] . These are reviewed 
below. 



A. Structure factors 

For lattice Hamiltonians, wavevectors that differ by 
the reciprocal lattice vectors G are equivalent. For the 
Fourier components of the local density of states p q this 
implies that Pq+c e (e) = p^ ttice (e) for any G. To un- 
derstand why this equivalence is not observed in experi- 
ments we must take into account the real space structure 
of the Wannier wavefunctions of electrons in the conduc- 
tion band. Here we study the effects of a single-band 
tight-binding model, in which Bloch states at the Fermi 
level can be written as a superposition of localized atomic 
orbitals Vk(r) = ^ R e !kR ^(r - R). 

We begin by projecting the Bogoliubov-de Gennes 
(BdG) wavefunctions in terms of the single-band wave- 
functions, 



u n(r) = ^a£Vk(r), 

k 



If we know how an operator © acts on the Bloch wave- 
functions Vk, ©V'k = X)k' Qkk'Vv, then the above rela- 
tion induces an action on through Oak = Y]^/ Ok'Ok'k 
(and similarly for &k). Thus, the solutions to the BdG 
equation 



will be independent of the Wannier wavefunction 0(r) 
once we determine the action of the BdG operator on 
the Bloch wavefunctions -0k- 

For positive biases e > (the e < case can be ana- 
lyzed analogously), the LDOS is given by 

< ys W= / dW^<(rK(r)<5(e-i;„) 

= 5(e ~ E n ) < n <,J(K k',q), 

n kk' 

where 

J(k,k',q) = J d 2 re^k(r)Vv(r) 

= ^r<5(q-(k-k') + G)x 
G 

^V k ' R j d 2 re 4qr 0*(r)(/)(r-R). 

R 

If we assume that the relevant electronic wavefunction is 
well localized, we can ignore terms involving the overlap 
across different sites (R ^ 0) in the last integral. Then, 
the only dependence of J on k and k' is through the 
crystal momentum conservation condition, and we find 



with 



C+V A+W 



< yS (e)-^qPq ttiCO ( £ ) ( 14 ) 

< ttice (e) = £%-£>ra£ +q 

k,n 

5 q = | d 2 r|0(r)| 2 e^ r . 



One immediately recognizes that in (14) p^ ttlcc (e) is the 
Fourier component of the lattice density of states that we 
analyzed in the earlier sections, and S q is the structure 
factor determined by the atomic wavefunctions. Peaks in 
the STM spectra arise from /Oq ttlcc (e), whereas S q only 
provides additional wavevector dependence. Hence in our 
tight binding model we expect that wavevectors which 
differ only by reciprocal lattice vectors have peaks at the 
same energies, but with generally different intensities. 



B. Current Normalization Condition 

An additional subtlety of STM experiments in [1-5] is 
the space-dependent normalization used. It is natural to 
assume that the tunneling matrix elements do not change 
appreciably with energy over the energy range of interest. 
Thus, if z is the height of the STM tip above the sample, 
and r is its 2D coordinate along the plane of the sample 
surface, then the differential tunneling conductance g can 
be written as 
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g(r,z,e) = f(r,z)p^ s (r,e), 

where p phys (r, e) is the 2D density of states in the CuO 
plane. The experiments in [1-5] adjust the z coordinate 
at every point r along the surface, so as to keep the cur- 
rent at Vf fixed at a predetermined value If . The differ- 
ential conductance normalized in this fashion is 

ff moas (r,6)=/(r)^ s (r, £ ), 

where /(r) = /// J^ Vf de p phys (r, e). Let us now discuss 
some properties of g meas (r, e). 

The spatial variation in /(r) is dominated by the in- 
homogeneous quasiparticle weight within a unit cell. To 
see this, write 

p phys (r, e) = p per (r, e) + p TS B(r, e), 

where p pe r(r, e) is periodic with the lattice and is of order 
one, whereas pTSB(r, e) breaks lattice translational sym- 
metry and is of order Vo in our formalism. If we define 
/per(r) = /// f* Vf dep peT (r,e), then 

/(r) = / per (r)(l - H(r)) 

in terms of a TSB function H (r) of order Vo- 
lt is convenient to absorb / por into p by introducing 
a new function, p'(r,e) = / por (r) j o phys (r, e). Due to the 
symmetry properties of / por , p' is simply related to p uttlcc 
through a modified structure factor 



~mcas / ,\ 
.Vq+GVV 



P q (e) = S q pif tice (e) 



P er(r)|0(r)|-e 



2 iqr 



Expressing g mcas in terms of p' , 



3 (e) 



P '^ e) -j (0^- k ^ (£) ' (15) 



we see that 5 q cas (e) gets "direct" contributions from 
structure in the LDOS at wavevector q, as well as 
"shadow" contributions from structure in the LDOS 
at other wavevectors k, whenever -ff q _k is non-zero. 
Whereas (15) is an exact relation, it is useful to trun- 
cate it to order Vq by keeping in the second term only 
those contributions coming from the neighborhood of the 
reciprocal vectors k « G, 



3 W 



Pq(e)-« q pLto co W + 0(K, 2 ) (16) 

^ Hq-G S' G (17) 



In this approximation the shadow contribution to 
<7™ eas (e) factorizes into the space-dependent factor a q 
and the space-averaged density of states /ojflg 00 oc 
(g mcas (e)). From (16) and (17) we can verify an impor- 
tant property of the tunneling spectra 



3 W 



J q+G 

SL 



when G is a vector of the reciprocal lattice and q is not. 
Hence, we expect <7 q eas (e) and ff^+cK 6 ) to nave peaks at 
the same energies but in general with different overall 
amplitudes. 

An interesting question to ask is whether it is possible 
to analyze experimental data in a way that would allow 
to separate direct and shadow contributions to g q oas (e). 
Below we demonstrate that this is possible using an exact 
sum rule obeyed by p q ttICC (e). Regardless of the model 
used and the nature of the symmetry breaking pertur- 
bation, the sum over all frequencies of p q ttlcc (e) should 
be identically zero for all q different from the reciprocal 
lattice vectors G: 

rfe^ ttice (e) = -^(0|{c k+qCT ,cL}|0) 

~°° ko- 

^2^(27r) 2 «5( q -G). 

G 

In principle, this identity can be used to remove the 
shadow contribution in equation (15). In particular, for 
q ^ G, combining the approximate result (16) and our 
knowledge of p l ^ lcc (e) from experiment, we can fix a q 
by requiring the sum rule to be obeyed: 



p q M = .C as ( e ) 



ffL x x .< cas (^ (18) 

mea.fi I ^\\ 



r Vmax mcas 
_ / mcas/ \\ J-Fmax ^ 
^ ^ '' f^MAX ; mcas 

J-14lAX 



>(e))de 



Here Vmax should be chosen sufficiently large so that the 
ratio of the two integrals is close to its saturated value, 
yet it should be small enough that we are still justified in 
using a single band model and a local picture of electron 
tunneling. 

For completeness, we also list two other sum rules 
obeyed by tunneling spectra. By construction, at every 
wavevector q ^ the function g q cas (e) must satisfy the 
normalization condition 



peVf 

Jo 



deg™ as (e)=0. 



One can derive an independent sum rule if we restrict 
the class of symmetry breaking Hamiltonians to effective 
one-particle operators (10) [This includes all perturba- 
tions considered in this work.] Then the e-weighted av- 
erage of /9 q ttlce (e) will be, for q ^ G, 

deepf^(e) = -]T<0| { [c k+qCT , H] , } |0) 
2 x 

kcr 

For the basis functions discussed in Section II we find 
that only the Vk q describing simple CDW gives finite 
contributions after summing over k. Hence, 
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f 

J —I 



rfeep' attice (e) 



2U 



cdw 



It is important to point out that this sum rule will be 
spoiled by shadow contributions in (15), and is only of use 
if these have been previously removed, using for example 
procedure in equation (18). 

As a useful consistency check of our formalism, one can 
easily verify that expressions (7), (8), and (9) satisfy both 
sum rules. We emphasize, however, that although these 
expressions are only correct to linear order in perturba- 
tion strength, the sum rules are non-perturbative and 
therefore hold to all orders in perturbation theory. Fur- 
thermore, their validity is not affected by the introduc- 
tion of finite quasiparticle lifetimes as, for any normal- 
ized symmetric distribution g(e), J de (a + /3e)g(e— e ) = 
a + /3e . By contrast, the average of /9 q (e) weighted by 
any other power of e is sensitive to details of quasiparticle 
smearing. 

Unfortunately, these sum rules are of limited imme- 
diate use, since the bulk of the integration comes from 
large energies, whereas current experiments only probe 
a relatively narrow range of biases about the chemical 
potential. 



VI. PHOTOEMISSION 

Before concluding, we would like to propose a way of 
identifying weak charge ordering in photoemission ex- 
periments that could supplement current STM studies. 
A common signature of a strong charge ordering in the 
ARPES experiments is the presence of shadow bands: the 
electron spectral function at momentum k acquires an 
additional peak at the energy of the quasiparticle at mo- 
mentum k + Q . For weak charge order the shadow bands 
may be difficult to observe: when the energy difference 
between £4. and -Ek+Q is large, mixing between quasipar- 
ticles is negligible and the intensity of the shadow peaks 
is vanishingly small. Strong mixing only occurs when 
states k and k + Q are nearly degenerate, although in 
this case the two peaks are hard to distinguish since they 
are close in energy. Thus we expect to observe an in- 
crease of the apparent lincwidth of quasiparticles when 
the latter satisfy the degeneracy condition = i?k+Q 
and are strongly affected by the charge order. For exam- 
ple, in the case of the Bi2Sr2CaCu20s+5 band structure 
shown on Figure 3 we expect an anomalous increase in 
the apparent quasiparticle linewidth at the points A, A', 
B, and B' on the Fermi surface. 



VII. CONCLUSIONS 

To summarize, we considered the effects of weak trans- 
lational symmetry breaking on the d-wave superconduct- 
ing state of the cuprates. For systems with periodic 



charge order we derived an explicit formula for the en- 
ergy dependence of the Fourier component of the local 
density of states for several types of order, including sim- 
ple charge density wave, electron kinetic energy and su- 
perconducting gap modulations. We argued that within 
a one band model the STM spectra observed in [1,2,5] 
cannot be explained by a simple charge density wave but 
require the existence of some form of (anomalous) dimer- 
ization, i.e. modulation in the electron hopping or in the 
superconducting pairing amplitude. We discussed a situ- 
ation in which charge order has finite correlation length 
due to pinning by impurities. In this case the LDOS 
has Fourier components for a range of momenta around 
the ordering wavevector Q. For different wavevectors q, 
peaks in p q (e) will occur at different energies, although 
the peak dispersion is weak, in agreement with [3,5]. We 
also considered systems in which translational symme- 
try breaking comes not from charge ordering but from 
impurities. We found that the Fourier components of 
the LDOS in this case have peaks for a wide range of 
wavevectors and strong dispersion of these peaks is con- 
sistent with the STM experiments of Refs. [3,4]. 
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FIG. 1. Order parameters with wavevector Q = (^f,0) 
considered in this paper. Figures (a) and (b) correspond to 
site and bond centered CDW respectively. Black circles corre- 
spond to sites of higher electron density, white circles to sites 
of lower electron density, and gray circles to sites with the 
average electron density. Figures (c) and (d) describe lon- 
gitudinal and transverse dimerizations respectively. Heavy 
lines correspond to bonds with higher tunneling amplitude, 
and dotted lines to bonds with lower tunneling amplitude. 
Anomalous dimerization may be shown schematically as on 
figures (c) and (d), with heavy and dotted bonds describing 
higher and lower pairing amplitudes. 
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FIG. 2. Energy dependence of the Fourier component of 
the local density of states pcj(e) at Q = (27r/4, 0) for vari- 
ous cases of charge ordering. Bi2Sr2CaCu208+5 type band 
structure is assumed. We show a direct calculation based on 
formulas (7,8,9). The curves correspond to the CDW (solid), 
longitudinal dimerization (dashed) and anomalous longitudi- 
nal dimerization (dotted) orders. To simplify the comparison, 
pq(e) is multiplied by —1 for CDW, and by i for anomalous 
longitudinal dimerization. In addition, subsequent curves are 
shifted vertically by 0.6. Results for both kinds of transverse 
dimerization are qualitatively similar to corresponding longi- 
tudinal results and are omitted for visual clarity. 
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FIG. 3. Fermi surface for Bi2Sr2CaCu208+,5 . Dashed lines 
correspond to the quasiparticles that satisfy -Ek+Q = £k and 
are strongly affected by ATI when Q = (27r/4, 0). Quasipar- 
ticles on line a hybridize with quasiparticles on line a' (and 
similarly for lines b and &', and curves c, c' and d, d'). Cross- 
ings of these lines with the Fermi surface (points A and A', B 
and B') give the minimal energy of such quasiparticles: 0.5Ao 
and 0.7Ao respectively. Contributions from these points pro- 
duce sharp peaks at energies 20 and 30 meV in Fig. 2. The 
Van Hove singularity for the Bogoliubov quasiparticles at en- 
ergy Aq leads to a peak at 40 meV. 
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FIG. 5. pg(e) for Bi2Sr2CaCu20s+,5 type band struc- 
ture, ordering wavevector Q = 27r/4 and a combina- 
tion of charge orders: longitudinal dimerization and CDW, 
1.05(CDW) + (long. dim.) (solid line); anomalous longitudi- 
nal dimerization and CDW, (anom. long. dim.) + 0.2(CDW) 
(long-dashed line). The same smearing is assumed as in Fig. 4. 
For clarity, curves have been offset vertically by 0.25. 
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FIG. 4. The results in Fig. 2 are shown after smearing over 
an energy range of 8 meV. To simplify the comparison, pq(e) 
was multiplied by —1 for CDW and anomalous transverse 
dimerization, and by a factor of — \ for transverse dimeriza- 
tion. 
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FIG. 6. YBa2Cu3 06.35 type band structure and ordering 
wavevector Q = (2tt/8, 0). The inset shows Pq(e) for CDW 
(solid line) , dimerization (dashed line) , and anomalous dimer- 
ization (dotted line) separately. Main figure has the lin- 
ear combination (aid) + 0.2(CDW), which was displayed for 
Bi2Sr2CaCu20s+5 in Fig. 5 (the other linear combination is 
nearly zero and thus omitted). The same smearing is assumed 
as in Fig. 4. 
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FIG. 7. Dispersion in the (0, 0) to (ir, 0) direction in a sys- 
tem with charge order with randomness (momentum is mea- 
sured in units of 2ir). Charge order is assumed to have Gaus- 
sian distribution centered around wavevector (2n/4, 0) with 
the width 2n/20. The function V(q) in (11) is shown, up to a 
scale, as the thick solid curve. For visual clarity, only results 
corresponding to the linear combination 1.05(CDW) + (ld) are 
displayed. Each curve corresponds to a different bias; starting 
from the bottom, the biases are 8 mV, 12 mV, 16 mV, and 
20 mV. Throughout, the quasiparticle smearing is fixed at 8 
meV. 
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FIG. 9. (a) Diagram of Hartree type used in computing 
the RPA response for a system of interacting electrons in the 
vicinity of a CDW instability. An external field (e.g. an im- 
purity potential) pins the CDW. (b) and (c) Contributions 
beyond the Hartree approximation, which may be subdomi- 
nant depending on the model used. Their inclusion leads to 
a frequency-dependent self energy. 
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FIG. 8. Dispersion in the (0, 0) to (n, 0) direction in a 
model with impurity induced quasiparticle scattering (mo- 
mentum is measured in units of 2n). Each curve corresponds 
to a different bias; starting from the bottom, the biases are 8 
mV, 12 mV, 16 mV, and 20 mV. Unlike other computations in 
this paper, the quasiparticle smearing is fixed at 2 meV. This 
is done since the main features in these curves are averaged 
out for the usual smearing of 8 meV. 
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